Transitions To the Long-Resident State in coupled chaotic oscillators 
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The behaviors of coupled chaotic oscillators before complete synchronization were investigated. 
We report three phenomena: (1) The emergence of long-time residence of trajectories besides one 
of the saddle foci; (2) The tendency that orbits of the two oscillators get close becomes faster with 
increasing the coupling strength; (3) The diffusion of two oscillator's phase difference is first enhanced 
and then suppressed. There are exact correspondences among these phenomena. The mechanism 
of these correspondences is explored. These phenomena uncover the route to synchronization of 
coupled chaotic oscillators. 

PACS numbers: 05.45.Xt. 



Synchronization is a universal and fundamental be- 
havior occuring in various fields pi Q. Differ- 
ent synchronous states can be observed in coupled 
chaotic systems such as complete synchronization (CS) 0> 
generalized synchronization(GS)0, phase synchro- 
nization (PS) (at! j measure synchronization (MS)[6j and 
so on. What people care about in studies of synchro- 
nization is often the critical coupling, in which the syn- 
chronization can achieve [Toj|. and the behavior near this 
threshold. The dynamic in regimes far from the global 
synchronous state is related to the mechanism that makes 
the synchronization be a stable state. In this paper we 
report three phenomena which can help us to understand 
the problem above. We also uncover a mechanism of syn- 
chronization in systems which have two or more saddle 
foci. 

When two or more chaotic systems were linearly cou- 
pled into a network, their differential equation can be 

represented by X = F(X) + eT <g> CX, where X = 



77 A< 



), and N is the node number of network. 



e denotes the coupling strength. T : W l — > K m char- 
acterizes the coupling schemes among the variables of 
the nodes in a network. C = M — D, where M de- 
notes the adjacency matrix of the network (the ele- 
ment Mij denotes the number of the edges that link 
node i and j), D is a diagonal matrix, and satisfies 
Da — Sj=i In this paper, we take the TV Lorenz 

systems (x = a(y — x),y = rx — y — xz, z = xy — f3z, 
where a = 10, r = 28, = 8/3) as our nodes and couple 
them into an array (the periodic boundary condition is 
applied) with Tij = = 1, 2, 3) except Tn = 1. 

There are three fixed points in phase space of a sin- 
gle Lorenz o scillator. Two of them are saddle foci 
(±y/[3(r — l),± Ak //3(r — l),r — 1) and another is a saddle 
node(0,0,0). The trajectory usually revolves around the 
two saddle foci alternatively. Fig^a) is the time serial 
of x-component of a single Lorenz system. The revolu- 
tion about one of the saddle foci can be regarded as one 
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FIG. 1: (a):Time series Xi(t) for a single Lorenz system. 
(b):Residence time statistics P(T). (c)(d):Coupled Lorenz 
system(iV = 2),e = 0.75. 



dynamic phase. Fig^b) is the statistics P(T) of the resi- 
dence time T that the system resides in a phase. This dis- 
tribution consists of several almost discrete peaks. When 
the two systems are coupled together, new phenomenon 
can be observed. FigQJc) is one of the time series of 
x-component in a system consisting of two Lorenz os- 
cillators. It can be observed that the residence time in 
certain phase is longer than the situation in Fig^a) . Ac- 
cordingly, a new wide continuous smooth peak appears 
in the distribution of residence time in Fig^d). The 
above phenomenon suggests a new state that we call long- 
residence state(LRS) in coupled chaotic systems. While 
the new peak in the distribution indicate the emergence 
of LRS, one can define P = ^ P(T)dT as the proba- 
bility of resident time that is larger than To. For a larger 
T , To is a measure of the emergence of LRS. FigCJa) is 
the relationship between e and -Po(T> = 10). Obviously, 
for a weak coupling, To ~ 0; But when the coupling e 
is greater than a critical point e\ « 0.36, To increase 
sharply, suggesting the emergence of the LRS. 
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FIG. 2: (a):P ~ e, To = 10. (b):Power spectrum of 
a;(t) for a single Lorenz system. (c):A(/) of xi(t) for coupled 
Lorenz system(iV = 2), e = 0.75. (d):f max ~ e, /o = 1. 
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FIG. 3: The four largest LEs for N = 2 Lorenz oscillators. 
There are two critical coupling strength e\ = 0.36 and £2 = 
1.38. 



On the other hand, the emergence of LRS can also 
be observed through the power spectrum of the chaotic 
motion. Fig[3fb) is the spectrum amplitude of a sin- 
gle Lorenz system, and Fig|2Ic) gives that of a coupled 
system with two oscillators. A new peak appears in the 
high-frequency regime. The relationship between this fre- 
quency peak in the region [/o, 00) and coupling strength 
e is plotted in FigEtd) ■ A critical coupling strength £1 
can also be found, indicating the emergence of LRS, and 
the critical point is same as that in Fig[2Ja). 

To uncover what happens near the critical point of 
LRS, it is instructive to study the Lyapunov expo- 
nent (LE) spectrum is useful. Figgis the four largest 
LEs of two coupled Lorenz oscillators. At £i=0.36 and 
£2=1-38, one of the zero exponents is found to become 
negative respectively. Therefore, the LRS is accompa- 
nied by a topological transition of the chaotic attrac- 
tor and is an intrinsic bifurcation embedded in com- 
plicated motion. One of the Lyapunov exponents be- 
coming negative means the decrease of the dimension 



of the chaotic attractor. On the other hand, when CS 
achieve at £2, the state of the coupled system is on an 
invariant synchronous sub-manifold, and the dimension 
of the attractor will become one half of that for uncou- 
pled system. Are there any relations between LRS and 
the change of the Lyapunov exponents. To answer the 
question, let's study the evolution of the trajectory dis- 
tance between two oscillators. We define the Euclidean 
distance as r(t) = y/ (x 2 - xi) 2 + (y 2 - y\) 2 + {z 2 - zi) 2 
. When the synchronization is achieved, r(t) will tend 
to rapidly. In the non-synchronous range, the behavior 
of r(t) is complex, and we prefer to studying its statisti- 
cal characters. By defining p(r)dr as the probability of 
the distance r located in r — > r + dr, the accumulative 
distribution P(R) = J Q p(r)dr describes the portion of 
trajectory distance that is smaller than R. In the non- 
synchronous regime P(R) — > for R — ► 0, and P(R) — > 1 
when R — > R m ax, where R m ax is the determined by the 
attractor size. When CS is achieved, one has P(R) = 1. 
FigEIa) shows the behavior of P(R) for different coupling 
strengths. It is interesting that P(R) obeys a power law 
for R < 1, i.e. P(R) oc R a . This indicates p(r) oc r"" 1 
, implying p(r) may change from a Gaussian-like func- 
tion to a Poisson-like function when a decreases below 
l(p(r) — > when r — > R ma x)- The variation of the ex- 
ponent versus the coupling e is shown in FigQJb). It is 
very interesting to notice that at £ < £1, a > 1, and 
when £ > £1, a < 1. This manifests the emergence of the 
long-resident state. 

the emergence of LRS implies longer rotations be- 
fore switching to another rotating saddle focus. This 
should be closely related to the phase dynamics of cou- 
pled Lorenz oscillators. In Lorenz system, a dynam- 
ical phase 9i = tg^^vLit) / 'u % x [t)\ can be introduced, 
where u l x (t) = Zi — (r — 1) and u l y (t) — \J x\ + yf — 
^2(3(r - l),(i = 1,2,- •• ,N). The definition of PS is 
that the average frequency of two or more oscillator is 
equal to each other, so the synchronization is in terms 
of the "average" of the phase. In fact, the phase defined 
above is not a well-defined quantity(no better definition 
of a single phase been proposed for oscillations with mul- 
tiple rotation centers). Therefore diffusion process can be 
observed for the phase difference between tow oscillators, 
i.e.AO = 6i — 9j,i,j = 1,2, • • • , N are time-dependent. To 
understand the changes of diffusive process, we calculate 
the second center moment of the difference (Ad 2 ) — (A9) 2 , 
where (■) means the average of ensembles. FigQJc) de- 
scribes the evolution of the difference at e = 0.75. It is 
approximately proportional to time, therefore this pro- 
cess is the Brownian motion. We can further define the 
diffusion coefficient, 

D =to* ^[(A6 2 ) - (A8) 2 } 

1 — >oo AL 

In Fig^Jd), we give the behavior of the diffusion coeffi- 
cient D against the coupling e. It can be found that the 
diffusion is first enhanced for weak couplings and then 
depressed for larger couplings. There is a peak near the 
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FIG. 4: AT = 2 (a):P(R) ~ i?. (b):a ~ e. (c):The behavior 
of {A6> 2 } - (A6») 2 . (d) The diffusion coefficient D against the 
coupling strengths. 
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FIG. 5: (a)(b):The spatiotemporal behavior of a lattice of 
N = 100, e = 0[(a)], e = 100[(b)]. (c)(d): The evolution of 
x(t) for two coupled Duffing oscillators( x + 0.5i + x 3 — x — 
0Acos(t) = 0) with different e. 



coupling £i where the LRS emerges. This correspon- 
dence implies that the appearance of LRS does favor to 
PS. Because there are two rotating centers in a single 
Lorenz system, a small couplin g wi ll do harm to the syn- 
chronization of coupled system[12j. When the system is 
in the LRS, oscillators usually prefer to staying in one of 
the phase, and the trajectories of two oscillators are easy 
to close. 

In fact, a system with more oscillators may exhibits 
stronger effect of LRS. Figs|5fa),(b) gives the spatiotem- 
poral patterns at e — and e = 100 respectively, for 
N = 100. The oscillator is labeled when it is in one 
scroll. One can find that the spatiotemporal pattern at 
e = is a random pattern, but obviously a very long stay 
of oscillators in one scroll can be found for e = 100. 

Finally, we want to stress that the LRS we observed 
here is generic. It can be observed in different situations. 
We have checked coupled Duffing oscillators(as shown in 
FigsJSfc) and (d), the residence time of two coupled Duff- 
ing oscillators is obviously longer than that of a single 
oscillator), Chua circuits, and Chen's oscillators All 
these systems exhibit the LRS behaviors for moderate 
couplings. Even if there are some parameter mismatches, 
a LRS can still be observed. Furthermore, we also find 
that the LRS does not depend on the coupling forms, 
e.g., x—,y— or z— coupling, and local or global. 

In conclusion, we explored the dynamics of coupled 
chaotic oscillators with multiple scrolls (saddle foci) prior 
to global synchronization. We find a generic transition 
to the long-resident state (LRS), where the oscillations 
experience a rather long duration in one scroll. This tran- 
sition is manifested by several different behaviors. First, 
one zero LE in LE spectrum become negative. Second, 
the transition to LRS is accompanied by the qualitative 
changes in the trajectory-distance distribution and power 
spectra. Third, we show that the emergence of the LRS is 
accompanied by the enhancement-depression transition 
of the diffusion of the phase difference. The phenomenon 
of LRS is generic, i.e., it is irrelevant to the dynamics, 
number, and the interaction types of oscillators. 
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